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1.  INTRODUCTION 


Recently,  a  useful  generalized  form  of  paraxial  wave  optics  has  been 

developed  that  car.  treat  paraxial  wave  propagation  through  any  optical  system 

1  2 

that  can  be  characterized  by  a  complex  ABCD  ray  matrix.  ’  It  has  been  shown 
how  Huygens'  integral  for  propagation  through  a  cascaded  series  of 
conventional  elements,  including  Gaussian-shaped  limiting  apertures  located 
arbitrarily  along  the  optical  Dath,  can  be  accomplished  in  one  step,  using  the 
appropriate  (complex)  ABCD  matrix  elements  for  that  system.  This  approach  has 
been  useful  for  handling  complicated  multielement  optical  resonators,  as  well 
as  resonators  with  the  useful  variable-reflectivity  mirrors  that  are  now  being 
developed. 1 

In  Ref.  2,  this  approach  has  been  generalized  to  the  case  of  propagation 
of  partially  coherent  light  and  link  propagation,  including  the  effects  of 
tilt  and  random  jitter  of  the  optical  elements  and  distributed  random 
inhomogeneities  along  the  optical  path  (e.g.,  clear  air  turbulence  and 
aerosols).  However,  the  results  presented  in  Ref.  2  for  propagation  in  the 
presence  of  random  inhomogeneities  were  limited  to  a  real  ABCD  system  (i.e., 
where  all  the  ray  matrix  elements  are  real).  Therefore,  the  results  of  Ref.  2 
are  limited  to  propagation  through  random  media  where  no  significant  apertures 
or  stops  are  located  between  the  input  and  output  planes  of  the  optical 
system.  Additionally,  for  propagation  between  conjugate  image  planes,  the 
overall  B  matrix  element  is  identically  zero  (for  real  ABCD  systems),  and  some 
of  the  expressions  in  Section  7  and  Appendix  C  of  Ref.  2  must  be  applied 
carefully  because  of  the  apparent  singularities  that  result.  We  will  show 
that  if  the  diffractive  effects  of  the  optical  elements  are  included  in  the 
analysis,  no  such  singularities  occur,  and  all  the  expressions  yield  finite 
results. 

Thus,  the  major  thrust  of  this  report  is  to  consider  propagation  through 
a  random  medium  in  a  general  complex  ABCD  paraxial  optical  system.  Explicit 
i  Ull3  ar-H  >)-r  i  vpc)  for  tn®  rorj.’  ting  second -order  statistical  moments  of 
the  complex  optical  field.  In  particular,  we  derive  a  general  expression  for 
the  mean  irradiance  of  a  point  source,  the  mutual  coherence  function,  the  log- 
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amplitude,  the  phase  and  wave  structure  function,  and  the  variance  of 
irradiance  for  weak  scintillation  conditions.  In  Section  2,  we  uiscuss  some 
general  features  of  the  scattered  field  through  terms  of  second  order  in  the 
fluctuating  part  of  the  index  of  refraction.  In  Section  3,  we  discuss  the 
propagation  of  a  point  source  and  obtain  an  explicit  expression  for  the  mean 
irradiance  of  a  point  source  propagating  through  a  complex  ABCD  system.  In 
Section  4,  we  present  explicit  results  for  various  second-order  moments  of  the 
complex  field.  In  Section  5,  we  present  some  specific  results  for  the  mean 
square  irradiance  of  a  spherical  wave  for  the  Kolmogorov  spectrum.  We  also 
give  explicit  results  for  the  conjugate  image  plane  and  Fourier  transform 
plane  propagation  geometries.  Finally,  in  Section  6,  we  consider  adaptive 
optics  through  real  ABCD  systems  and  derive  an  expression  for  the  mean 
irradiance  of  an  adaptive-optics  laser  transmitter  through  such  systems. 

The  analysis  in  this  report  and  in  Ref.  2  applies  primarily  to  optical 
systems  that  contain  apertures  of  Gaussian-shaped  variations  in  amplitude 
transmission  across  the  optic  axis.  However,  the  analysis  also  applies,  at 
least  as  a  first  approximation,  to  any  transversely  varying  system  whose 
transmission  has  a  quadratic  variation  to  first  order  near  the  optic  axis.1 
Thus,  any  "soft"  aperture  whose  transmission  varies  with  transverse 
coordinate,  at  least  near  the  axis,  in  the  approximate  form 

T(x)  =  T0(  1  -  a2x2/2) 


where  TQ  is  the  on-axis  transmission  and  the  coefficient  *s  given  by 

.  j  (d2T(x) 

2  “  "  T  1  .  2 

o  dx 

can  be  approximated  to  first  order  by  a  Gaussian  aperture  and  thus  by  a 

complex  ABCD  matrix.  As  a  result,  the  complex  paraxial  analysis  is  a  good 

approximation  as  long  as  the  resulting  beam  wave  remains  sufficiently  close  to 

2 

the  axis  so  that  |a ^x  /2 |  <<  1  across  the  main  portion  of  the  beam. 
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In  contrast  to  propagation  through  real  ABCD  systems,  the  statistical 
moments  of  the  field  that  result  after  propagating  through  a  complex  ABCD 
system  are  not  spatially  homogeneous.  That  is,  the  statistical  moments  depend 
explicitly  on  the  variables  of  interest  rather  than  on  the  corresponding 
difference  variables.  This  dependence  occurs  because  the  introduction  of 
finite  limiting  apertures  destroys  the  spatial  symmetry  of  the  field  of  a 
point  source.  For  example,  after  propagating  through  a  complex  ABCD  system, 
the  mean  irraaiance  distribution  of  a  point  source  is  no  longer  uniformly 
distributed  over  a  sphere;  rather,  it  is  limited  in  space  (i.e.,  it  becomes  a 
beam  wave).  As  a  result,  all  statistical  moments  of  the  complex  field  are,  in 
general,  explicit  functions  of  absolute  position  variables  rather  than  of 
corresponding  difference  variables. 

In  this  report,  we  consider  the  case  where  the  intervening  space  between 
optical  elements  is  nearly  that  of  free  space  w.th  a  refractive  index  given  by 
n  =  1  +  n^.  We  assume  the  random  quantity  to  have  zero  mean  and  a  root 

mean  square  value  much  less  than  unity.  The  present  analysis  is  thus 
restricted  to  a  weakly  inhomogeneous  medium  with  an  electrical  conductivity 
equal  to  zero  and  a  magnetic  permeability  equal  to  one.  We  further  assume 
that  the  characteristic  scale  length  of  the  inhomogeneities  is  much  greater 
than  the  optical  wave  length  and  that  the  characteristics  of  the  medium  do  not 
change  appreciably  during  a  period  of  oscillation  of  the  electromagnetic 
field.  For  example,  at  near  IR  and  optical  wavelengths,  this  condition  is 
satisfied  in  the  atmosphere.  Additionally,  we  consider  the  propagation  of 
scalar  monochromatic  fields  and  omit  the  explicit  time  dependence  in  the 
formulation  presented  below. 

Finally,  for  simplicity  in  notation,  we  consider  cyl indr i cal ly  symmetric 
optical  systems,  because  the  extension  to  orthogonal  systems  is  straight¬ 
forward.  In  cyl  indr  ical ly  symnetric  optical  systems,  the  ray  transfer  matrix 
M  is  a  2  »  2  matrix  given  by 
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where  A,  B,  C,  and  D  are  complex  numbers.  For  the  case  considered  here,  where 
the  input  and  output  planes  are  in  Free  space 


Det  (M)  =  AD  -  BC  =  1 


(1.2) 
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2.  GENERAL  CONSIDERATIONS  FOR  THE  SCATTERED  FIELD 


In  the  intervening  space  between  the  optical  elements,  the  scalar  wave 
equation  is 


V2U  +  k2n2(r)U  =  0 


(2.1) 


where  U  is  a  typical  component  of  the  field,  k  is  the  optical  wave  number,  and 


n(r)  :  1  +  n ^  ( r ) 


(2.2) 


where  <n^>  =  o  and  |n^|  <<1.  We  indicate  the  statistical  average  by  angular 
brackets  and  the  spatial  coordinate  r  =  (r,z),  where  z  is  assumed  to  be  along 
the  optical  axis  and  r  is  a  two-dimensional  vector  transverse  to  the  optical 
axis.  As  indicated  in  Fig.  1,  we  assume,  without  loss  of  generality,  that  the 
input  plane  is  located  at  z  :  0  and  the  output  plane  is  at  the  plane  z  =  L. 

We  employ  the  Rytov  transformation,^  which  consists  of  setting  U  =  exp  [  tp  ] 
in  Eq.  (2.1)  and  thus  obtaining 


2 

V  Hi 


;  V4-)  +  k  [  1 


n  i  ( r ) 


]  --  0 


It  is  customary  to  seek  a  solution  for  ip  as  a  power  series  in  e1  =  (2n1  +  n^ ) 


*  ~~  l.  K 


m  =  0 


where  ipQ  is  zero  order  in  e  1 ,  ^  is  first  order  in  and  so  forth. 

In  this  report,  we  include  terms  through  second  order  in  e ^ .  The  second- 
order  Rytov  approximation 


U(r)  =  exp[ipo(r)  +  ik  1  ( r )  +  ip^ir)] 


(2.3) 
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INPUT  PLANE  OUTPUT  PLANE 


Fig.  1.  Schematic  Representation  of  Propagation  through  an  Arbitrary  ABCD 
Optical  System 
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which  contains  all  terms  through  second  order  in  e  ^ ,  is  the  lowest  order  non¬ 
trivial  approximation  that  conserves  energy  to  the  order  of  the  approximation 

and  -hat  gives  the  correct  average  field  and  the  correct  phase  and  intensity 

a 

statistics. 

Now  UQ  =  exp  4/  is  the  field  in  the  absence  of  the  inhomogeneities  and  is 

p 

given  In  the  paraxial  approximation  by 


Uq(£o,L)  =  -  exp  ( -ikL)  J  d2r  U.(r)  exp  [-  ||  (Dp2  -  2r-£  +  Ar2)j  (2.4) 

where  is  the  field  in  the  initial  plane  and  A,  B,  and  D  are  the  ray  matrix 

elements  of  the  overall  system  between  the  input  and  output  planes.  The 
quantity  L  appearing  on  the  right-hand  side  of  Eq.  (2.4)  is  the  optical  path 
length  of  a  ray  traveling  along  the  optic  axis.  Because  ol  the  approximate  ns 
used,  we  assume  that  this  quantity  is  equal  to  the  total  propagation  distance 
between  the  input  and  output  planes. 

To  within  an  arbitrary  constant,  the  general  solution  for  4^  is  given  by 


0 


where  U1  is  the  first-order  Born  approximation 


U^r)  =  -k‘ 


d3r. 


G(r, r. 


e1(rl)Uo(rl) 


(2.5) 


(2.6) 


where 


e 1 ( r )  =  2n  j ( r )  +  n2(r)  (2.7) 

The  integration  in  Eq.  (2.6)  extends  over  the  region  of  space  where  n^(r)  is 
different  from  zero,  and  G(r,r1>  is  the  appropriate  Green  function  in  the 
absence  of  the  inhomogeneities. 

For  direct  propagation  between  transverse  planes  located  at  z .  and  Zg, 
the  Green  function,  in  the  paraxial  approximation,  is  given  by 
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"  ■  1  "  1  >  >  ■  'I  I  ■  I  I  ]  II  II  ■  II  I  II.  ,1 .  LB  I  I . 

f 

i 

|  G(r1.?2)  =  G(r1  -  ?2) 

2 

exp  [-ik(z_  -  z  )]  ik(r  -  r  ) 

-  -  -7-7 - g-,-1-  exp |  -  -77^ - =V| 

4h(z2  -  z^  F1  2(z2  - 

®  In  the  presence  of  a  complex  paraxial  ABCD  system,  this  Green  function  is  now 

given  by  9 


G( r ! ’ r2)  =  - 


exp 


-ik(z2  -  z 1 ) 


4nB 


exp[  - 


ik 


Z1Z2 


2B 


(A  r.  -  2r- r»  +  D  r  )  j 
z,z_  1  - 2  z,z_  2  1 

Z1Z2  12  12  (2.8) 


where  we  have  employed  the  notation  A  to  represent  the  A-matrix  element 

z  ]  >  z2 

for  propagation  between  z ,  and  z?  and  similarly  for  the  quantities  B„  and 

\  C  Z.  y  Z*2 

D_  _  . 

z,z2 

The  second-order  Rytov  approximation  is  given  by^ 


Up(r) 

*P(r)  =  — -  -  5*:(r)  (2.9) 

UQ(r)  d  1 

where  4'1  is  the  first-order  Rytov  approximation.  The  quantity  U2,  the  second- 
order  Born  approximation,  is  given  by 


U2(r)  =  (-k2)2Jd3r1Jd^r2G(r,r1)e1(r1)G(r1,r2)e1(r2)  UQ(r2)  (2.10) 

where,  for  the  complex  paraxial  system,  UQ  is  given  by  Eq.  (2.4)  and  G  is 
given  by  Eq.  (2.8) . 

Equation  (2.3)  can  be  used  to  compute  all  second-order  statistical 
moments  of  the  field,  as  obtained  in  a  plane  transverse  to  the  optic  axis  at 
propagation  distance  L,  through  terms  of  second  order  in  the  fluctuations.  We 
can  show  that  all  such  second  moments  can  be  obtained  from  three  basic  field 
moments . 
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r 

t 


F-|  =  <4s2(g)  > 

F2  r 

and 

f3  =  < 4--,  (£1 )  <1*!  (E2)  > 

where 


o 


(2.11a) 

(2.11b) 


(2.11c) 


(2.12) 


and  Hi.  is  given  by  Eq.  (2.5). 

For  example,  the  average  field  <U;  =  UQ  exp(<4>  >),  the  mutual  coherence 
function,  is  given  by 


r(£r£2)  -  <U(g  T  )U  (£2)> 


“  exp[  v  <t>2 ( £ i )  > 


* 

<4>2(£2)>  + 


(2.13) 


and  the  normalized  variance  of  irradiance  (under  weak  scintillation 
conditions)  is  given  by 


<I2(£)> 

<I(£)>2 


-  exp  [ 4x2(£) ]  -  1 


where 


x2(£)  =  |  Re(F2  ♦  F3) 


(2.14) 


and  Re  denotes  the  real  part. 
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To  obtain  the  field  in  the  output  plane  from  an  arbitrary  field  in  the 
input  plane  (e.g.,  via  a  Huygens-Fresnel  integral),  it  is  sufficient  to 
calculate  the  quantities  F^,  F2,  and  F^  for  the  special  case  where  the  sources 
of  the  fields  are  point  sources.  To  calculate  these  quantities,  we  denote, 
by  ^(jD.r)  and  ^(fiiE),  the  first-  and  second-order  Rytov  approximation  to  the 
optical  field  at  (g,L)  due  to  a  point  source  located  at  (r,0),  respectively. 
The  calculation  of  the  statistical  quantities  F ^ ,  F2,  and  F^  is  given  in  the 
Appendix,  with  the  results  that 

F1  ;  c«2(£,r)> 

=  — itk2  J  dz  Jd2K®  (K;z )  (2.15) 

0 

F2  =  <<'1(er£1)l'i(e2’-2)> 

=  2nk2  J  dz  J  d2K  ®  (K  ;2)exp(-K2e. /k) 

0  n  1 

x  exp  (-iK-{[g  Re(B0z/B)  +  £  Re  (B^/B)] 

♦  2i  [P  Im(B0z/B)  ♦  R  Im(BzL/B))})  (2.16) 

f3  =  ^1(21,r1)®1(£2,r2)> 

=  -  2nk2  J  dz  J  d2K  ®  (K;z)  exp(iK2e/k) 

0  n 

x  exp  {- iK -  f e( b0z/b)  +  £(BzL/B)H  (2.17) 

where  ®n(K;z)  is  the  three-dimensional  spectrum  of  the  index  of  refractive 
fluctuations  at  propagation  distance  z  evaluated  at  Kz  =  0.  The  quantities 
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Bq2,  BzL,  and  B  e  BqL  are  the  B-matrix  elements  for  propagation  through  the 
system  from  0  to  z,  z  to  L,  and  0  co  L,  respectively 


8 


BA  b  , 
Oz  zL 


(2.18) 


8  =  lm(e),  where  Im  denotes  the  imaginary  part,  and 


£  =  -  fi2  (2’ 19a) 

P  =  \  (fi1  +  22)  (2- 19b> 


r  =  r 1  -  r2  (2. 19c) 

and 

-  =  5  (-1  +  -2)  (2. 1 9d ) 

In  general,  Bq2  and  BzL  (and  hence  8  )  depend  on  z,  depending  on  the  optical 
system  under  consideration.  Through  these  quantities,  the  optical  system 
affects  the  second-order  statistical  moments  of  point  sources.  For  the 
important  case  of  an  isotropic  power  spectrum  (e.g.,  the  Kolmogorov  spectrum), 
d2K  >  2^KdK  and 


exp  (-iK-V)  -  JQ(K | V | ) 

where  Jq  is  the  zero-order  Bessel  function  of  the  first  kind  and  V  is  the 
(complex)  vector  within  the  braces  in  Eqs.  (2.16)  and  (2.17).  In  the  rest  of 
this  report,  we  will  be  dealing  with  spatially  isotropic  power  spectra. 

Examination  of  Eq .  (2.15)  reveals  that,  through  terms  of  second  order  in 
the  fluctuations,  the  quantity  F1  =  <$2>  =  <U2>/UQ  is  real  and  is  independent 
of  the  observation  point  and  the  matrix  elements  of  the  optical  system.  In 
particular,  Eq.  (2.15)  is  identical  to  what  is  obtained  for  line-of-sight 
propagation  over  a  distance  L.  As  indicated  in  Eqs.  (2.16)  and  (2.17),  the 
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quantities  F2  and  F^  depend  on  the  specific  properties  of  the  optical  system 
and  are  spatially  nonhomogeneous ,  because  they  depend  explicitly  on  the 
variables  r  ^ ,  £2,  ,  and  £>2*  However,  for  the  special  case  of  real  ABCD 

systems,  F2  and  F^  are  spatially  inhomogeneous,  because  there  are  no  limiting 
apertures  that  inhibit  the  propagation  of  spherical  waves  through  the 
medium.  Note  that  for  1  ine-of-sight  propagation,  Bqz/B  =  z/L,  Bz^/B  = 

(L  -  z),  and  8  =  z(L  -  z)/L.  In  this  case,  Eqs.  (2.16)  and  (2.17)  become 
identical  to  the  corresponding  results  given  in  the  literature. ^ 

Thus,  Eqs.  (2.15)— (2. 17)  contain  all  the  effects  of  the  arbitrary  complex 
paraxial  ABCD  system  on  the  second-order  statistical  properties  of  point 
sources.  They  are  now  the  basis  for  calculating  all  the  second-order 
statistical  moments  of  interest. 


14 


3.  MEAN  FIELD  AND  IRRADI ANCE  OF  A  POINT  SOURCE 


Consider  a  point  source  located  at  the  origin  of  coordinates,  which 
radiates  a  total  power  W.  To  within  a  multiplicative  phase  factor,  we  can 
show  from  Eq .  (2.4)  that  the  field  in  the  observation  plane  in  the  absence  of 
the  inhomogeneities  is  given  by 


U  (p)  ^  g)  exp  (-ikL  - 

°  U,|Bf 


(3-D 


from  which  it  follows  that  the  irradiance  IQ  =  |UQ|  is  given  by 


I  (p)  =  -  —5  exp  [-kj  Im(D/B)  jp2  ]  (3.2) 

°  |  Bp 

p 

For  real  ABCD  systems,  1  =  W/4*B  ,  which,  for  1 ine-of-sight  propagation, 

2 

reduces  to  I  =  W/4tiL  . 

Examination  of  Eq.  (3.2)  reveals  that  in  the  presence  of  a  complex 

paraxial  ABCD  optical  system,  the  irradiance  distribution  of  a  point  source 

2 

has  a  Gaussian  shape  with  a  1/e  spot  radius  pQ  given  by 


po  =  k| Im(D/B) I 


(3.3) 


As  an  elementary  example,  consider  a  point  source  located  at  distance  to 
the  left  of  a  Gaussian  lens  of  real  focal  length  f  and  1/e  transmission 
radius  o  (see  Fig.  2).  If  the  observation  plane  is  in  the  right-hand  focal 
plane  of  the  lens,  the  results  are 


2z.  i 

B  =  F  (1  -  -V) 


z.  2z.i 

D  =  1  _  —  _  — — 
f  ,2 
ko 


(3.4) 


(3.5) 
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r 


INPUT  PLANE 


LENS  OUTPUT  PLANE 


POINT 

SOURCE 


Fig.  2. 


Propagation  Geometry  of  a  Point  Source  Located  a  Distance 
Left  of  a  "Thin  Gaussian"  Lens  of  Real  Focal  Length  f  arid 
Transmission  Radius  a.  The  output  plane  is  assumed  to  be 
right-hand  focal  plane  of  the  lens. 


z1  to  the 
1/e2 
in  the 
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If  we  substitute  Eqs.  (3.4)  and  (3.5)  into  Eq.  (3-3),  the  result  is 

p2  =  (fo/z.)2  +  (2f/ko)2  (3.6) 

o  1 

p 

The  square  of  the  1/e  beam  radius  is  a  sum  of  two  terms:  the  first  term,  on 
the  right-hand  side  of  Eq.  (3.6),  is  a  geometric  magnification;  the  second 
term  gives  the  diffraction  effects  of  the  finite  lens. 

When  inhomogeneities  are  present,  the  mean  field  is  given  by 

<.U(p)>  =  Uo<exp[ii»1(g)  +  ^(R))'  (3.7) 

If  g  =  ln(y)  is  a  random  variable  through  terms  of  second  order  in  the 

p 

fluctuations  with  mean  <g>40  and  mean  square  <g  >,  then  correct  through  terms 
2 

of  the  order  of  e ^ , 

<y>  =  exp  [<g>  +  |  <(g  -  <g>)2>]  (3.8) 

This  result  may  be  used  to  obtain  the  average  of  the  right-hand  side  of  Eq. 

2  2 

(3-7).  Hence,  through  terms  of  the  order  of  e^  (»  4n^), 

<U(p)>  =  Uq  exp  ( < 4>2 > )  (3-9) 

where  UQ  is  given  by  Eq.  (3.1)  and  <<{>2  >  is  given  by  Eq.  (2.15).  Apart  from 
the  UQ  term,  the  average  field  of  a  point  source  is  identical  to  that  obtained 
for  1 ine-of-sight  propagation  through  the  random  medium. 

The  mean  irradiance  of  a  point  source  is  obtained  as 

<I(r)>  =  <jU(p)|2>  (3.10) 

If  we  substitute  U(p)  =  UQ(g)  expti^p)  +  ^(r)]  into  Eq.  (3.10)  and  apply  Eq. 
(3-8)  to  obtain  the  average,  the  result  is 

<  1  (r)>  =  IQ(R)  exp  [  Mg)  ]  (3.11) 
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where  IQ(p)  is  given  by  Eq.  (3.2)  and 

A(£)  =  2Re<4>2>  ♦  <|^1(e)|2>  (3.12) 

If  we  substitute  Eqs.  (2.15)  and  (2.16),  with  r^  =  £2  =  0  and  2i  -  E2  =  £» 
into  Eq.  (3.12),  we  obtain  the  following  equation  for  an  isotropic  index  of 
refraction  power  spectrum 

L  as 

A(p)  =  4*2k2  f  dz  /  dK  K*  (K;z) 

0  0  n 

[exp(-K26./k)I0(2o.pK)  -  1]  (3.13) 

where 

a.  =  Im(B0z/B)  (3.14) 

and,  in  Eq.  (3.13),  IQ(x)  is  the  zeroth-order  Bessel  function  of  the  first 
kind  of  imaginary  argument.  Equations  (3.10)  and  (3.13)  are  the  general 
expressions  for  the  mean  irradiance  profile  of  a  point  source  propagating 
through  an  aribtrary  complex  ABCD  system  in  the  presence  of  random  media,  as 
specified  by  *n(K;z).  For  real  ABCD  systems,  A(p)  =  0,  as  it  should  because 
of  spatial  synmetry.1*  Because  Ig(x)  2  1,  the  contribution  to  A(p)  from  the 
Bessel  function  term  is  always  greater  than  zero,  indicating  that  the 
inhomogeneities  broaden  the  irradiance  pattern. 

For  small  values  of  p,  one  can  readily  obtain  a  power  series  expansion  of 
A(p)  as 

L  CD 

A ( p )  =  -4*2 k2  /  dz  J  dK  K*  (K;z) 

0  0  n 

x  [1  -  exp(-K26i/k)  -  a2p2K2  exp  (-K^/k)*. . .  ]  (3.15) 
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Consider  the  Kolmogorov  spectrum  in  the  inertial  subrange 


3 


n 


(K;z) 


=  0.033  C2(  z)  K“11/3 


(3.16) 


2 

where  C  (z)  is  the  index  structure  constant  profile.  If  we  substitute 
n 

Eq.  (3.16)  into  Eq.  (3.15)  and  perform  the  integration  over  spatial  frequency, 
the  result  is 


A<p)  -  -4.25  k 


7/6  J  dz  c^zHe^z)]576 


+  3.65  k2p2  J  dz  C2(z)  a2(z)  [k/B.(z)]1/6-. ..  (3.17) 

0 

2 

For  the  propagation  geometry  illustrated  in  Fig.  2,  with  C  ^  0  between  the 

o  n 

point  source  and  the  lens  only,  and  7-  >>  k  o“,  we  readily  obtain 


A(p)  -  -  ( 1.62o/po)573  [1  -  1.25  (ko/2f)2p2--..)  (3.18) 


where 


pq  =  [1.46  k 


0 


dz  C2(z)(z/Zl)5/3]'375 


(3.19) 


is  the  spherical  wave  lateral  coherence  length  of  a  spherical  wave  for 
propagation  between  z  =  0  and  z  =  Zy  From  Eqs.  (3.11)  and  3.18),  we  see  that 
the  first  and  second  term  on  the  right-hand  side  of  Eq.  (3.18)  give  the 
inhomogeneities-induced  reduction  in  the  on-axis  mean  irradiance  and  beam 
broadening . 
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4.  GENERALIZED  SECOND-ORDER  MOMENTS 


In  this  section,  we  present  general  expressions  for  various  second-order 
statistical  moments  that  arise  in  propagation  through  complex  ABCD  systems  in 
a  random  medium. 

4.1  CORRELATION  FUNCTIONS 

We  have 


*,(£,£)  =  x(£.r)  -  iS(g, r)  (4.1) 

where  x(Ri£)  is  the  log-amplitude  and  S(g,r)  is  the  phase  at  a  point  2  in  the 
output  plane  resulting  from  a  point  source  at  r  in  the  input  plane.  For 
simplicity  in  notation,  we  omit  subscripts  on  both  x  and  S.  We  can  then  show 
that 


® ^ ^ E.  1 » P-2  ’  — 1’— 2^  =  x  x(.P^  ) x(£2  2^  x  r  2  Re(F2  +  F3^  (4.2) 

bs(£1.22’  -1’-2)  =  <S(Et2i)S(£2’-2):>  =  5  Re(F2  '  F3J  (4-3) 

and 

BSx(Er£2;  -1’-2)  =  <s(e1»-1)x(E2’-2)>  =  5  In,(F2F3)  (4.4) 

where  B  is  the  log-amplitude,  B^  is  the  phase,  and  B„  is  the  phase  log- 
X  ox 

amplitude  correlation  functions. 

Except  for  adaptive  optics,  discussed  in  Section  6,  the  general  two-point 
to  two-point  moments  are  not  of  general  interest,  although  they  can  be 
obtained  directly  by  substituting  Eqs.  (2.16)  and  (2.17)  into  Eqs.  (4.2)- 
(4.4).  We  restrict  our  attention  to  the  special  case  of  a  single  point  source 
located  in  the  input  plane  at  r  =  0.  Following  the  above  procedure,  we  find 


21 


B  _(p.  ,p_)  =  2n2k2  /  dz  /  dK  K*  (K;z)exp(-KV/k) 

^  0  0  n 

k2b 

x  {Re[JQ(K  jgar+2iPa  .  |  )  ]  T  cos  [—— )  JQ(Kpar)}  (4.5) 


where  the  t  sign  refers  to  the  log-amplitude  and  phase  correlation  functions, 
respectively,  and 


v  ■  Re<>WB> 


B_  =  ReB 
r 


(4.6a) 

(4.6b) 


Similarly 


Be  (ElE-,)  =  2^k2  I  dz  /  dK  K$  (K;z)  expt-K^B. /k) 
“x  0  0 


k2b. 


{lm[J0(K|ear  +  2iPa.|)]  ♦  sinf-j^)  JQ(Kpor)}  (4.7) 


For  the  special  case  where  =  -g2»  Eqs.  (4.5)  and  (4.7)  become 

L  co 

B  _(p,-p)  =  2it2k2  J  dz  J  dK  K*(K;z)exp(-K26./k) 
X  i  x*  0  0  n  1 

K26. 


[l  ?  cos(-~ -))J0(2Kpar) 


(4.8) 


and 


B0  (p,-p)  =  2«2k2  f  dz  f  dK  K$  (K;z)  exp(-K2B./k) 
S*  0  0  n  1 


K2B 


sin(  — )  JQ(2Kpar) 


(4.9) 


The  variance  of  log-amplitude,  phase,  and  phase  log-amplitude  is  obtained  from 
Eqs.  (4.5)  and  (4.6)  by  setting  p,  =  p2-  ^inci 
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<X2>| *  Bx  S(£’£) 

<s  >'  x’ 


=  2n2k2  f  dz  f  dK  K®  (K;z)  exp( -K^B . /k) 
0  0  n  1 


L  L 


K2B 


I0(2Kpa.)  i  cos( — ip) 1 


(4.10) 


and 


<SX>  =  Bs  (fi,g) 


r  2*2k2  /  dz  J  dK  K®  (K;z)  exp(-K^0./k) 

0  0  n  1 


sin  (- 


k2b 


(4.11) 


4.2  MUTUAL  COHERENCE  FUNCTION 

The  mutual  coherence  function  (MCF)  is  defined  as 


r(£r£2)  =  ^(fil)U  (£2): 


(4.12) 


We  omit  the  deterministic  contribution,  UQ  (g^  UQ  (£2),  to  the  M(-F- 
contribution  can  be  obtained  directly  from  Eq.  (3.1).  We  now  have 


This 


U(£)  =  expO^g)  +  i|>2(£)  ] 


(4.13) 


where  41.  is  given  by  Eq.  (2.5)  and  is  given  by  Eq.  (2.9).  If  we  substitute 

I  d 

Eq.  (4.13)  into  Eq.  (4.12)  and  use  Eq.  (3.8)  to  obtain  the  statistical  average 

of  the  right-hand  side  of  Eq .  (4.12),  the  result  is 


# 


r’(£1,£2)  =  exp[2<®2>  +  <4'1(£1)4'1(E2)>] 


(4.14) 


I 


I 


! 
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If  we  .substitute  Eqs.  (2.15)  and  (2.16)  into  Eq.  (4.14),  we  obtain 

L  00 

r(£.,£p)  =  exp(-4*2k2  J  dz  J  dK  K$  (K;z) 

1  *  0  0 

*  [1  -  exp(-K2d./k)  J0(K|£ar  -  2Pa . | ] }  (4.15) 

In  general,  r  is  complex  and  spatially  inhomogeneous  with  respect  to  £1  and 
£2-  Equation  (4.15)  is  the  general  expression  for  the  MCF  of  a  point  source 
located  at  the  origin.  The  equation  applies  for  a  general  complex  ABCD 
paraxial  system,  as  specified  by  a  and  8.  For  real  ABCD  systems,  Eq.  (4.15) 
becomes  identical  to  the  corresponding  results  of  Ref.  2. 

For  some  applications,  knowledge  of  the  complex  degree  of  coherence  is 
required.  This  quantity  is  given  by  the  normalized  MCF^ 

>(E1.e2)  = 

If  we  substitute  Eq.  (4.14)  into  Eq.  (4.16),  we  obtain 

y(£1,£2)  =  exp{*,j(£1)>n2(£2)>  -  |  t  <  1 4*-,  (£1 )  |2>  +  <  l*>1(£2)  |2>1  >  (4.17) 

2  L  ®  P 

=  exp{-2irk  J  dz  J  dK  K*  (K;z)  exp(-K  8./k) 

0  0  n  1 

x  [ IQ(2a1P1K)  ♦  I0(2a.p2K)  -  2J0(Klpar  +  2iPa.|)]}  (4.18) 


r(£r£2) 

[r(E1,E1)rrE2’£2)1 


(4.16) 
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5.  I'  '(ADI ANCE  STATISTICS 


2 

For  weak  scintillation  (i.e.,  x  1).  it  is  well  known  chat  the 
irraaiance  statistics  can  be  expressed  as  second-order  moments  of  the 
field. ^  Here  we  present  a  general  expression  of  the  normalized  variance  of 
irradiance  of  a  point  source  located  at  the  origin 

2  I2<£>' 

°7  -  5  -  1  (5.1) 

1(2) 

where 

1(2)  =  | U( 2> |2  (5.2) 

and 

U(£)  =  Uo(£)  exp[  i1  j ( p )  +  4>2(2)]  (5.3) 

If  we  substitute  Eqs .  (5.2)  and  (5.3)  into  Eq.  (5.1),  making  use  of  Eqs.  (4.1) 
and  (3-8),  we  obtain  through  terms  through  second  order  in  the  fluctuations 

Oj(fi)  -  exp[4- ^(2) '-I  -  1  *  4<x2(2)>  (5.4) 

where  x^p)  is  given  by  Eq .  (4.10).  Effects  that  cause  beam  broadening  in 

the  output  plaiiu  (e.g.,  wave  front  tilt)  are  contained  identically  in  <I^>  and 
2  2 

<I'  .  Thus,  the  normalized  variance  is  independent  of  such  effects  and 
reflects  the  focusing/defocusing  properties  of  the  random  inhomogeneities  only 
along  the  optical  path. 

For  example,  consider  a  point  source  located  at  distance  Zj  to  the  left 
of  a  lens  of  real  focal  length  f  and  Gaussian  radius  o.  We  will  determine  the 
on-axis  irradiance  variance  at  a  distance  z?  to  the  right  of  the  lens.  We 
will  assume  that  the  intervening  space  between  the  lens  and  point  source  is 
filled  with  a  turbulent  medium,  which  is  characterized  by  the  Kolmogorov 
spectrum . 
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For  p  =  0  and  the  Kolmogorov  spectrum  [Eq.  (4.10)],  we  obtain 


I 


where 


<X2(0)  =  0.652  k2  /  dz  C2(z)  f  dK  K"8/3  exp(-K26./k) 
0  n  0  1 

x  [1  -  cos(K2B  /k) ] 
r 

=  0.326  k7/6  /  dz  C2(z)  b  (z)5/6W(z) 

0  n  r 


(5.5) 


(5.6) 


W(z)  =  J*  -yy J7  (1  -  cos  x)  exp[-y(x)x] 

0  x 

=  r(-5/6) -  (1  +  u2)8/12  cos[5/6  arctan( 1/u)  ] )  (5-7) 

B^z) 

u< z >  '  (Tu)  '5.8) 

and  r(  . .  . )  is  the  Gamma  function.^  For  real  ABCD  systems,  u  =  0  and  Eq.  (5.6) 
yields  the  results  given  in  Ref.  2.  For  y  <<  1,  Eq.  (5.7)  yields 
W(z)  -  -  r(-5/6)cos(5n/12)  =  1.73,  and  for  u  >>  1,  Eq.  (5.7)  yields 
W(z)  -  -  (5/ 12) r C -5/6)  u~7/6  -  2.78  y"7/6. 

Equations  (5.6)-(5.8)  are  the  general  expressions  for  the  on-axis  log- 
amplitude  variance  for  the  Kolmogorov  spectrum.  These  equations  are  valid  for 
an  arbitrary  complex  ABCD  system,  as  specified  by  6(z).  The  quantity  W(z)  is 
an  additional  weighting  function  that  results  from  the  optical  system.  In 
contrast  to  the  results  for  real  ABCD  systems,2  no  singularities  result  for 
complex  ABCD  systems. 

To  be  specific,  we  consider  two  well-known  examples  of  practical 
interest:  imaging  and  Fourier  transform  propagation  geometries. 


» 
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5.1  IMAGE  PLANE  VARIANCE  OR  I RRADI ANCE 
Here  we  have 


1_  J_  _  1 

Z1  *  Z2  "  f 

2 

For  simplicity,  we  assume  that  C  *  0  only  between  the  point  source  and  the 

2 

lens.  It  is  straightforward  to  show  that 


B  =  - 


2iZ1Z2 

ko2 


(5.9a) 


and  for  z  <  z 


EL  -  z 
Oz 


(5.9b) 


and 


zz„  2iz,z0 

bzl  •  ^  -  -H  <’  -  n> 


(5.9c) 


Thus,  for  z  <  z 


"r  *  Z  (’  -  t  > 


(5.  10a) 


and 


2,  2 

.  z  k° 

8i  •  2 


2z 


(5.10b) 


Hence,  from  Eq .  (5.8),  we  obtain 


a(z/z^ ) 

i  -  i 


(5.11) 


where 


a  = 


ko 
2z . 


(5.12) 
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Over  the  range  of  integration  in  Eq.  (5.6),  where  the  integrand  is  non-zero, 

z  -  z-|,  and,  hence,  y  ~  a.  Thus,  when  the  point  source  is  in  the  far  field  of 
2 

the  lens  (z1  >>k  o  ),  y  >>  1  over  the  range  of  integration.  Conversely,  for 

2 

near-field  conditions  (z1  <<,  ko  ),  y  >>  1  over  the  range  of  integration.  As  a 

result,  the  asymptotic  limits  for  the  on-axis  log-amplitude  variance,  for 
2 

constant  C^,  are  given  by 

<x2>  =  0.564  k7/6  C2  J  1  dz  6  (z)5/6 
n  0  r 


0.  124  k7/6  C2  z]1/6,  z  >>  ko2 
n  1  ’  1 


(5.13) 


and 


<x2>  =  0.906  k7/6  C2  J  1  dz  6  (z)5/6(y(z)  )“7/6 
n  0  r 


0.612  k7/6  cV^^^z^ko2)7^,  z1  <<  ko2 


(5.14) 


Equation  (5.13)  is  just  the  spherical  wave  variance  obtained  for  1  ir.e-of-sight 
propagation  over  a  distance  Zy,  whereas  Eq.  (5.14)  is  proportional  to  this 
variance  multiplied  by  an  averaging  aperture  factor  (z^k  o2)7/^  <<  1.  We  can 
explain  these  results  by  noting  that,  for  weak  scintillation  conditions,  the 
transverse  irradiance  correlation  length  for  propagation  over  a  distance  Zy  is 
of  the  order  (z-j/k)  .  For  an  imaging  geometry,  light  incident  on  the  lens 

from  the  left  converges  toward  the  image  point  on  the  right.  Then,  for  far- 
field  conditions  [o  <<  (z^/k)^],  the  entire  lens  lies  within  a  single 
correlation  patch,  and  no  aperture  averaging  is  expected.  For  near-field 

1  /p 

conditions  [o  >~>  (z^/k)  ],  many  independent  correlation  patches  are 

contained  over  the  imaging  lens,  resulting  in  a  reduced  spherical  wave 
variance  at  the  image  point.  For  the  Kolmogorov  spectrum,  this  reduction  or 
aperture  averaging  factor  is  proportional  to  the  -7/3  power  of  the  number  of 
independent  correlation  patches  over  the  aperture. ^ 
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2 

Equation  (5.6)  has  been  integrated  numerically  for  arbitrary  a  =  ko  /2z.. 

2  2  ' 
In  Fig.  3,  we  present  a  plot  of  <x  >  /<x  >  „  as  a  function  of  a. 

H  cl~U 

5.2  FOURIER  TRANSFORM  PLANE  VARIANCE  OF  IRRADIftNCE 

Here,  we  consider  the  case  where  the  point  source  is  in  the  left-hand 
focal  plane  and  the  observation  point  is  in  the  right-hand  focal  plane  of  a 

lens  of  real  focal  length  f  and  Gaussian  radius  o.  It  is  straightforward  to 

2 

show  that 


and  for  z  <  f 


B  =  f 


(1  - 

ko 


(5. 15a) 


B 


Oz 


z 


(5.15b) 


and 


BzL  8  B 


2zfi 

ko2 


Thus,  for  z  <  f,  we  obtain 


(5- 15c) 


B 


r 


,<2 

l  d 


+ 


(i 


1  +  a 


2 


(5. 16a) 


and 


where 


8i 


2 

az 


f(1  «■  a2) 


(5.16b) 


a  =  ko  /2f 


Hence,  from  Eq.  (5.8),  we  obtain 


u  = 


a(z/f ) 


2  , .  z. 

a  +  ( 1  -  ^) 


(5.17) 
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Normalized  Variance  of  Irradiance  as  a  Function  of  "a"  for  Image 
Plane  and  Fourier  Transform  Plane  Propagation  Geometries .  The 
quantity  a  =  ka  /2z  for  the  image  plane  geometry  and  ko  /2f  for  the 
Fourier  transform  plane  geometry. 


Examination  of  Eq.  (5.17)  reveals  that  in  both  the  near  field  and  far  field  of 
the  transform  lens,  p  =  6i/6r  <<  1.  In  the  far  field,  6r  =  z(1  -  |) ,  whereas 

in  the  near  Held,  b  -  z.  As  a  result,  the  asymptotic  limits  lor  the  Fourier 

r  2 
plane  on-axis  log-amplitude  variance  are  given  by  (for  constant  C^) 

2  n  nh  i  7/6  r2  11/6  ,2  /ci  a  \ 

>  =  0.124  k  C^z^  ,  z1  >>  ka  (5.18) 

and 

<x2>  =  0.307  k7/6  C2z]1/6,  z1  <<  ko2  (5.19) 

Equations  (5.18)  and  (5.19)  reveal  that,  in  the  far  field  (near  field)  of 

the  transform  lens,  the  on-axis  variance  of  log-amplitude  is  given  by  the 

spherical  wave  (plane  wave)  variance  for  direct  propagation  ever  a  distance 

z i . ^  In  contrast  to  the  image-plane  log -amplitude  variance,  no  aperture 

averaging  effects  are  obtained  in  the  transform  plane  for  near-field 

conditions.  These  asymptotic  results  can  be  explained  as  follows.  In  the  far 

field,  the  wave  impinging  on  the  lens  is  nearly  planar,  hence,  the  resulting 

field  on  the  right-hand  side  of  the  lens  is  approximately  a  spherical  wave 

converging  towards  the  focal  point.  Furthermore,  in  the  far  field,  a  << 

1  /? 

(z/k)  ,  and  no  aperture  averaging  is  expected.  Conversely,  in  the  near 

field  of  the  transform  lens,  the  wave  on  the  right-hand  side  of  the  lens  is 
nearly  collimated  (i.e.,  a  plane  wave  propagating  parallel  to  the  optic 
axis).  As  a  result,  the  light  that  arrives  on  the  optic  axis  in  the  transform 
plane  results  from  a  small  region,  of  characteristic  size,  of  the  order 
(f/k)  ,  about  the  axis  of  the  lens.  That  is,  the  light  arriving  at  a  point 

on  the  optic  axis  in  the  transform  plane  is  nearly  planar  and  emanates  from  a 
region  on  the  order  of  a  single  correlation  patch  in  the  plane  of  the  lens. 
Thus,  no  aperture  averaging  effects  are  to  be  expected. 

2 

For  arbitrary  values  of  a  =  ko  /2f,  Eq.  (5.6)  has  been  integrated 
numerically.  The  results  are  plotted  in  Fig.  3.  These  results  indicate  mild 
aperture  averaging  for  a  ~  1. 
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6.  ADAPTIVE  OPTICS  AND  ABCD  SYSTEMS 


In  this  section,  we  present  some  results  concerning  the  performance  of 
adaptive-optics  systems  in  an  ABCD  paraxial  optical  system  between  the  input 
and  output  planes.  For  simplicity,  we  consider  real  ABCD  systems,  although 
the  methods  outlined  here  can  be  used  directly  to  obtain  the  corresponding 
results  for  complex  ABCD  systems.  To  successfully  implement  adaptive-optics 
correction  to  either  a  laser  transmitter  system  or  to  an  imaging  system, 
accurate  measurements  must  be  made  of  the  phase  errors  associated  with  the 

7 

appropriate  propagation  path.  Then,  according  to  the  principle  of 

reciprocity,  if  the  appropriate  correction  (i.e.,  the  negative  of  the  measured 

errors)  is  applied,  the  performance  of  the  system  will  improve 

significantly.  We  consider  conventional  phase-conjugate  adaptive-optics 

systems,  where  a  beacon  signal  provides  an  estimate  of  the  desired  phase 

information.  We  assume  that  the  beacon  signal  is  a  spherical  wave  that 

originates  from  some  point  in  the  output  plane.  Sometimes,  the  beacon 

location  originates  at  the  object  of  interest  (i.e.,  the  aim  point  of  a  laser 

transmitter  or  the  object  of  the  imaging  system),  and  significant  adaptive- 

optics  improvement  can  be  expected.  Conversely,  for  some  applications,  the 

object  of  interest  is  moving  relatively  fast.  Therefore,  if  light  emitted 

from  the  object  is  used  as  a  beacon,  the  adaptive-optics  improvement  is 

limited,  because  the  beacon  signal  will  propagate  through  a  portion  of  the 

atmosphere  different  from  that  where  adaptive-optics  correction  is  desired. ^ 

These  and  other  effects  result  in  the  various  anisoplanatic  degradation 

7  fi 

discussed  in  the  literature.  ’ ° 

To  be  specific,  we  consider  an  adaptive-optics  laser  transmitter  and  seek 
to  determine  the  resulting  mean  irradiance  obtained  at  some  point  £i  in  the 
output  plane,  where  the  beacon  phase  information  originates  from  a  point 
source  located  at  some  other  point  £2  the  output  plane.  For  real  ABCD 
systems,  the  mean  irradiance  at  £^  is  given  by^ 


<I(£1)>  = 


1 2^B  1 2  /  ^  exP("  rc(£) 


(6.1) 
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where 


K(r)  =  /  d2R  U.(R  +  ^r)U*(R  -  |r)  exp(-  ^  r-R)  (6.2) 

# 

rc(r)  =  <exp[*c(r1,E1;B2)  +  «i»c(r2,£1  ;b2)  ] >  (6.3) 

IL  is  the  input  wave  function,  A  and  B  are  the  corresponding  matrix  elements 
for  propagation  through  the  entire  ABCD  system,  and 

*c(£tE;E2)  =  “  iSB(E.E2)  (6.4) 

is  the  "corrected"  turbulence- induced  wave  function.  In  Eq.  (6.4) 

4>(r,g)  =  'P1  ( r.g)  +  it»2(r,£)  (6.5) 

is,  through  terms  of  second  order  in  the  fluctuations,  the  Rytov  approximation 
to  the  field  at  r  in  the  input  plane  resulting  from  a  point  source  at  g^  in 
the  output  plane,  and  Sg(r,g2)  is  the  phase  of  the  beacon  field  at  r  resulting 
from  a  point  source  located  at  g2  in  the  output  plane.  The  laser  and  beacon 
are  assumed  to  have  the  same  wavelength.  Note  that  in  phase-conjugation 
adaptive-optics  imaging  systems,  part  of  the  optical  transfer  function  imposed 
on  the  system  resulting  from  the  random  medium  is  also  given  by  r  . 

For  a  given  input  wave  function  and  ABCD  system,  the  deterministic 
quantity  K,  given  by  Eq.  (6.2),  can  be  obtained  in  a  straightforward  manner. 

In  the  remainder  of  this  section,  we  concentrate  on  evaluating  the  statistical 
quantity  r  ,  the  "corrected"  mutual  coherence  function.  If  we  assume  that  the 
beacon  phase  is  a  2ero-mean  random  variable,  application  of  Eq.  (3.8)  to  Eq. 
(6.3)  yields  through  terms  of  second  order  in  the  fluctuations  that 

rc  =  exp  t “<  1 4<c |2 >  +  F(£rE2,E1  ;E2)  1  (6.6) 
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where 


F(£r£2’£r£2)  =  <4lc1(£r£1’E2)  "'c1(£2’fil;E2)> 

and 

^cl(£.e;E2)  =  ’‘’I ~  iSB<£>22)  (6-8) 

In  obtaining  Eq.  (6.6M6.8),  we  have  used  the  fact  that  for  real  ABCD  systems 

2Re<4>2>  +  (4^  |2  =  0 

which  can  be  obtained  directly  from  Eqs.  (2.15)  and  (2.16)  for  2  =  £  =  0.  We  have 
also  used  the  fact  that  from  statistical  stationarity ,  <x(£1 »Ei)SB(r1(22)>  = 

< X (£2 » £•]  )^g( £2 >E2)  ^  >  and  <S(  r^  ,2-j  )Sg(  r.|  ,22 )  ^  =  ^S( r2 , 2 )^g(£2*E2 )  • 

Upon  substituting  Eqs.  (6.7)  and  (6.8)  into  Eq.  (6.6),  neglecting  the 
correlation  between  phase  and  log-amplitude,  we  can  be  show  that 

rc  =  exp  [-^Dw(r)  -  |ds  (r)  +  d(r,2)l  (6.9) 

B 

where 

Dw<£)  =  <[x(£rE1)  ~  x(£2-E1))2>  +  <ts^£i»E1)  "  s(£2*E!^2>  (6.10) 

is  the  wave  structure  function  of  a  point  source  located  at  2i  in  the  output 
plane . 8 


(6.7) 


Ds  (r)  =  <(SB(rv22)  "  SB(r2,22)]2>  (6.11) 

B 

is  the  phase  structure  function  of  the  beacon  point  source  located  at  22  in 
the  output  plane,  and 
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d(r,£)  =  <[S(r1,£1)  -  S(r2,£1)j[SB(r1,£2)  -  Sg( r2 ,£2) ] : 
=  |  (Dss  (r,£)  +  (-£,£)  -  2D<,q  (£)] 


(6.12) 


where 


Dss  <£<£)  =  <|S(£1,£1)  -  S(r2,£2)]‘ 
B 


(6.13) 


is  the  structure  function  for  point  sources  located  at  £1  and  £2  in  the  output 
plane  evaluated  at  the  points  r^  and  £2  in  the  input  plane  [in  Eq.  (6.12), 
DSSg^  =  DSSb(0,P)].  The  quantity  d  is  the  correlation  function  between  the 


B 


residual  phase  difference  of  the  object  wave  across  the  input  plane  and  the 
beacon  wave  across  the  input  plane.  Because  they  are  statistically 
stationary,  Dw  and  Dg  are  independent  of  the  location  of  the  point  source,  in 
contrast  to  d,  which  is  a  function  of  £  and  £.  If  the  beacon  is  collocated 
with  the  object,  and  the  measured  beacon  phase  is  a  perfect  estimate  of  the 
object  phase,  rc  =  exp(-^  Dft),  where  is  the  log-amplitude  structure 
function.  This  results  in  the  best  possible  performance  that  can  be  obtained 
with  conventional  phase  conjugation  adaptive-optics. 

The  structure  functions  that  appear  in  Eqs.  (6. 10)— (6. 13)  can  be  obtained 
directly  from  the  results  of  Sections  2  and  <4.  By  substituting  Eqs.  (2.16) 
and  (2.17)  into  Eqs.  (4.2)  and  (4.3),  we  can  show  that  for  isotropic  power 
spectra 


D  (£)  =  8*2k2  J  dz  J  dK  K*  (K;z)  [l  -  J  (Ko.r)] 
0  0  n  U 


(6.14) 


D  (r)  =  4i2k2  J  dz  J  dK  K*  (K;z)  |l  -  J  (Ka  r) 


0  0 


I  1  ♦  cos(yj] 


O'  1 


(6.15) 
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and 


D-„  (r)  =  4n2k2  J  dz  J*  dK  K4  (K;z)[  1  +  cos(^p)  ] 


k2a> 


ss 


x  [1  -  J0(K-al£  + 


(6.16) 


where  r  =  r 


£  =  £i  "  H2»  E  =  Bl>£2 


a1  (z) 


BzL(2) 

B 


,  .  B0z<z> 
a2  2  ‘  B 


(6.17a) 
(6. 17b) 


and 


B(z)  = 


BOz(z)  BzL(2) 
B 


(6.17c) 


Equations  (6 .  14 )-( 6 . 17)  and  Eq.  (6.9)  represent  the  general  solution  for  the 

corrected  MCF  in  a  real  ABCD  system  for  use  in  calculating  the  adaptive-optics 

mean-irradiance  profile  when  correlations  between  phase  and  log-amplitude  are 

neglected  (the  usual  case  considered  in  the  literature).  If  this  latter 

restriction  is  relaxed,  it  can  be  shown  that  r  contains  an  additional 

c 

multiplicative  factor  of 

C0SlBS  -  Bs  <r2,ri;El,£2)l 


where  B  ,  the  correlation  function  of  the  beacon  phase  and  object  point  log- 

OgX 

amplitude,  is  given  by 


3S  x(^r£2;Er£2)  =  Z”2*2!  d z  J  dK  K*n(K;z) 


0  0 


sin  (i^)jo(K|air  +  a22|) 
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Within  the  framework  of  the  present  theory,  the  inclusion  of  the  above  term 
in  Tc  provides  the  coupling  between  phase-amplitude  and  log-amplitude  in 
adaptive-optics  systems.  For  certain  applications  (e.g.,  propagation 
conditions  where  the  interaction  between  atmospheric  turbulence  and  thermal 
blooming  is  important),  the  inclusion  of  this  term  is  necessary  to  obtain  a 
self-consistent  analysis. 

As  an  example,  we  consider  the  Kolmogorov  spectrum  in  the  inertial 
subrange^  for  propagation  conditions  where  geometrical  optics  is  valid.  The 
cosine  terms  appearing  on  the  right-hand  side  of  Eqs.  (6.15)  and  (6.16)  can  be 
replaced  by  unity.  As  a  result,  we  obtain 

rc(r)  =  exp[-Sc(r)]  (6.18) 

where 

S  (r)  =  2.91  k2  Jdz  C2(z)  ((a?r)5/3  ♦  (a.p)5/3 
0  ^ 

-  \  la2-  +  “^l573  -  |  l«2£  *  5/3l  (6.19) 

For  line-of-sight  propagation,  =  1  -  z/L  and  <*2  =  z/L,  which,  when 
substituted  into  Eq.  (6.19),  yields 

S  (r )  =  2.19  k2  J  dx  C2(z){[r(1  -  z/L) J5/3  ♦  (ez)5/3 
0 

-  ~  | r ( 1  -  z/L)  +  ezj5/3  -  |  | r( 1  -  z/L)  -  0z|5/3}  (6.20) 

where  0  =  g/L  is  the  vector  angular  separation  between  the  beacon  and  the 
object  point.  Equation  6.20  is  identical  to  previous  results  found  in  the 
literature  for  couinon  or  angular  anisoplanatism. ® 


7.  CONCLUSIONS 


t 


Within  the  framework  of  the  paraxial  and  the  second-order  Rytov 
approximation,  we  have  derived  explicit  general  expressions  for  calculating 
all  second-order  statistical  moments  of  an  optical  wave  propagating  in  a 
random  media  through  an  arbitrary  complex  paraxial  ABCD  system.  These 
expressions  are  extensions  of  the  corresponding  results  that  apply  both  to 
line-of-sight  propagation-*  and  propagation  through  real  ABCD  systems.  The 
approach  in  this  report  (i.e.,  Huygens-Fresnel  formulation  of  the  field)  uses 
complex  ABCD  matrices.  This  approach  provides  the  most  general  form  for 
analyzing  propagation  of  any  generalized  paraxial  wave  through  a  complex 
paraxial  optical  system  in  the  presence  of  random  inhomogeneities  distributed 
arbitrarily  along  the  optical  path.  According  to  the  Huygen's  integral 
approach,  an  arbitrary  optical  wave  propagated  through  an  arbitrary  complex 
paraxial  system  in  a  random  medium,  including  all  diffraction  effects  (even 
those  resulting  from  the  "Gaussian"  limiting  apertures),  can  be  accomplished 
in  one  step,  with  knowledge  only  of  the  ABCD  matrix  elements  of  the  system. 
Even  though  Huygen's  integral  has  been  used  extensively  in  the  literature,  its 
general  range  of  validity  is  perhaps  not  as  well  understood  as  it  ought  to  be. 
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APPENDIX:  DERIVATION  OF  F , ,  F2,  AND  F3 

In  this  appendix,  we  outline  the  derivation  of  the  basic  statistical 
quantities  F^,  F2,  and  F3  introduced  in  Section  2.  The  method  used  is  similar 
to  that  given  in  Ref.  9.  The  basic  difference  is  that  the  zeroth-order 
paraxial  Green  function,  rather  than  the  free  space  function,  is  given  by  Eq. 
(2.8). 

When  the  optical  field  in  the  absence  of  the  inhomogeneities  is  point 
sources,  the  field  U  ,  as  obtained  from  Eq .  (2.4),  is  given  by 


C—  x  rw^. 

I6 

Uo(-1 ,z)  =  “bI -  6XP 

Oz 


(Dn  r.  -  2r  -r  ♦  An  r2)] 


2Ba  '“Oz  1 
Oz 


( A- 1 ) 


where  C1  =  -ik/2n  and  the  point  source  is  located  at  (r,0). 

A. 1  DERIVATION  OF  F1 

If  we  substitute  Eqs.  (2.8)  and  (A-1)  into  Eq.  (2.10),  taking  the 
statistical  average  and  rearranging  terms,  we  obtain 


C.k  -  ikL  ,  ^ 

<U2(E’L,>  *  -  --72-  !  d  r,  1  d S  "c'VV'VA.L1 


x  H(r  ,?  )  exp  [- 


(A0z2r  -  2-  -2  +  D0z2r2}l 


0z2  2 


x  exp  j-  ig—  (A  r  -  2r  •£  -  D  )] 
z}L  1  1 


( A-2) 


where 


Bc(r1  “  r2)  =  <E1(r1)E1(r2)> 


is  the  correlation  function  of  the  fluctuations  of  the  dielectric  function 


[note  that  Bn  =  <n(r^)n(r2)>  «  4B£ ]  and 
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=  <Bz,z /' 


ik 


2  1 


2B 


Z2Z1 


(A  r2  -  2r0 • r  +  D  r2 
1  Z2Z1  2  "2  '  2Z1 


2)J  ( A— 3 ) 


Next,  we  expand  in  terms  of  its  two-dimensional  Fourier  transform¬ 
er  ”  r2)  =  J  d2K  (K,z i  -  z2)  exp[-iK-(r1  -  £2)] 


(A-4) 


The  quantity  F£(K,az)  is  non-zero  only  for  |  Az  |  -  scale  length  of  the 
inhomogeneities.  We  assume  that  the  largest  scale  length  is  much  smaller  than 
any  propagation  distance  of  interest  (i.e.,  the  smallest  distance  between 
optical  elements).  Consider  the  integration  over  Z2.  We  only  get  a  non-zero 
contribution  to  this  integral  for  |z^  -  z^ |  ♦  0.  Thus,  Az^z^  =  Dz^z  =  1,  and 


z  ,z 


1^2 


).  Hence,  the  function  H  (r^,r2)  is  given  by 


H  = 


lim 

z2z1  -  0 


-ik(r.  -  r  ) 
[exp( - J - )] 


Z1Z2 


2ui  r,  . 

IT  s(ii  -  £z> 


(A-5) 


where  6(r)  is  the  two-dimensional  delta  function.  To  obtain  this  limit,  we 
assume  that  k  has  a  small  imaginary  part,  i.e.,  the  wave  propagates  in  a 
slightly  absorbing  medium.  Substituting  Eqs.  (A-4)  and  (A-5)  into  Eq.  (A-2) 
and  noting  that  the  z-dependence  of  the  exponential  factor  is  slowly  varying 
with  respect  to  that  of  F  ,  we  obtain 


:U2(e,L)>  = 


C1k4(2iri/k)e_lkL 
(4n ) 2 


‘  (b0z,bj,l)'' 


J  d2K  J  dz.  J  d2r . 
0  1  1 


ik  ( *  2  ^  ^  2' 

55—  (»0z,r  •  *  Vl- 

Oz  ^  1 


exp(-  5B~T  (Az,Lrt  -  2VP  +  Dz,LP'")1  J  dZ  Fe(^’z)  (A'6) 
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where  Z2  =  z-|  in  the  exponential  terms  on  the  left-hand  side  of  Eq.  (A-6), 
and,  because  F  -  0  for  z  2  scale  length  of  the  fluctuations,  the  upper  limit 
on  the  integral  over  z  can  be  extended  to  infinity.  From  Ref.  3,  we  have 


f  F  (K,2)dz  =  7i®  (K) 

0  e  "  E 

=  4*®n(K)  ( A-7) 


where  <J>n < K_ )  is  the  three-dimensional  spectral  density  of  the  index  of  refrac¬ 
tion  fluctuations  evaluated  at  Kr  =  0. 

Substituting  Eq.  (A-7)  into  Eq.  (A-6),  we  obtain 

<U?(£,L)>  ;  C  k4(2iTi/k)elkL  J  dz  J  d2K®  (K ) I  (A-8) 

0  n 

where 


and 
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Because  MqL  ,  M^M^,  where  Mz  is  the  ABCD  matrix  for  propagation  from  z^ 
to  z 2 ,  we  can  show  that  AzLBq2  +  bzld0z  =  B0L  and>  thus 

J  =  J  d  r1  exp  l-  y  -  lk  (§—  +  g— )  •r1  ] 

Oz  zL 

=  (2*i/W)  (^~^)  exp  [i§£  {—  +  ~)  ]  ( A  —  1 1 ) 

0L  Oz  zL 

where  B  is  given  by  Eq .  (2.18).  To  ensure  convergence  of  this  integral,  we 
have  again  assumed  that  Imk  <  0.  Combining  Eqs.  (A-10),  (A-11),  and  (A-8)  and 
using  the  relations 


a 
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we  finally  obtain 


CU_ (E,L ) >  =  -U ■  (£,L)  Ttk2  /  dz  /  d2K4> 
*  0 


(K;z) 


(A-12) 


where 


C  e 

Uo(£’L)  =  exP  I"  ^  (A0Lr2  '  +  D0Lp2)J  (A"13) 

is  the  field  of  the  point  source  in  the  absence  of  the  inhomogeneities,  and  we 
have  tacitly  assumed  a  slowly  varying  dependence  of  *n  on  propagation 
distance.  Thus,  from  Eqs.  (A-12),  (A-13),  (2.12),  and  (2.11a),  we  obtain  the 
results  given  by  Eq.  (2.15). 

A .  2  DERIVATION  OF  F2  AND  F-, 

The  derivation  of  both  F2  and  F^  follows  exactly  the  same  steps  that 
result  in  Eq.  (19)  of  Ref.  9,  except  that  Eq.  (2.8)  is  used  as  the  zeroth- 
order  Green  function  instead  of  the  free  space  Green  function.  The  final 
result  of  such  a  procedure  leads  directly  to  Eqs.  (2.16)  and  (2.17). 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  "architect-engineer”  for 
national  security  projects,  specializing  in  advanced  military  space  systems. 
Providing  research  support,  the  corporation's  Laboratory  Operations  conducts 
experimental  and  theoretical  investigations  that  focus  on  the  application  of 
scientific  and  technical  advances  to  such  systems.  Vital  to  the  success  of 
these  investigations  is  the  technical  staff's  wide-ranging  expertise  and  its 
ability  to  stay  current  with  new  developments.  This  expertise  is  enhanced  by 
a  research  program  aimed  at  dealing  with  the  many  problems  associated  with 
rapidly  evolving  space  systems.  Contributing  their  capabilities  to  the 
research  effort  are  these  individual  laboratories: 

Aerophyalcs  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
spacecraft  structural  mechanics,  contamination,  thermal  and  structural 
control;  high  temperature  thermomechanics,  gas  kinetics  and  radiation;  cw  and 
pulsed  chemical  and  exclmer  laser  development  including  chemical  kinetics, 
spectroscopy,  optical  resonators,  beam  control,  atmospheric  propagation,  laser 
effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  state-spectf ic  chemical  reactions  and 
radiative  signatures  of  missile  plumes,  sensor  out-of -f ield-of -view  rejection, 
applied  laser  spectroscopy,  laser  chemistry,  laser  optoelectronics,  solar  cell 
physics,  battery  electrochemistry,  space  vacuum  and  radiation  effects  on 
materials,  lubrication  and  surface  phenomena,  thermionic  emission,  photo¬ 
sensitive  materials  and  detectors,  atomic  frequency  standards,  and 
environmental  chemistry. 

Computer  Science  Laboratory:  Program  verification,  program  translation, 
performance-sensitive  system  design,  distributed  architectures  for  spaceborne 
computers,  fault-tolerant  computer  systems,  artificial  intelligence,  micro¬ 
electronics  applications,  communication  protocols,  and  computer  security. 

Electronics  Research  Laboratory:  Microelectronics,  solid-state  device 
physics,  compound  semiconductors ,  radiation  hardening;  electro-optics,  quantum 
electronics,  solid-state  lasers,  optical  propagation  and  communications; 
microwave  semiconductor  devices,  mlcrowave/mllllmeter  wave  measurements, 
diagnostics  and  radiometry,  microwave/millimeter  wave  thermionic  devices; 
atomic  time  and  frequency  standards;  antennas,  rf  systems,  electromagnetic 
propagation  phenomena,  space  communication  systems. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metals, 
alloys,  ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  non¬ 
destructive  evaluation,  component  failure  analysts  and  reliability;  fracture 
mechanics  and  stress  corrosion;  analysis  and  evaluation  of  materials  at 
cryogenic  and  elevated  temperatures  as  well  as  in  space  and  enemy-induced 
environments. 

Space  Sciences  Laboratory:  Magnetospheric,  auroral  and  cosmic  ray 
physics,  wave-particle  interactions,  magnetospheric  plasma  waves;  atmospheric 
and  ionospheric  physics,  density  and  composition  of  the  upper  atmosphere, 
remote  sensing  using  atmospheric  radiation;  solar  physics,  infrared  astronomy, 
infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  ionosphere  and  magnetosphere; 
effects  of  electromagnetic  and  particulate  radiations  on  space  systems;  space 
instrumentation. 


